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Abstract 

We present a determination of the QCD parameter A in the quenched 
approximation (nj = 0) and for two flavours (n/ = 2) of light dynamical 
quarks. The calculations are performed on the lattice using 0(a) improved 
Wilson fermions and include taking the continuum limit. We find A^^^'^q = 
259(l)(20)MeV and A^^2 = 261(17)(26) MeV, using tq = 0.467 fm to 
set the scale. Extrapolating our results to five fiavours, we obtain for 
the running coupling constant at the mass of the Z boson a"'^(m^) = 
0.112(1)(2). 
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1 Introduction 



The parameter A is one of the fundamental quantities of QCD. It sets the scale for 
the running coupling constant and it is the only parameter of the theory in 

the chiral limit. Usually A is defined by writing as{fi) as an expansion in inverse 
powers of ln(yU^/A^). For such a relationship to remain valid for all values of fi, 
A must change as flavour thresholds are crossed: A —>■ A„^, where Uf indicates 
the effective number of light (with respect to the scale fi) quarks. 

A lattice calculation of A requires an accurate determination of a reference 
scale, the introduction of an appropriate non-perturbatively defined coupling, 
which can be computed accurately on the lattice over a sufficiently wide range of 
energies, as well as a reliable extrapolation to the chiral and continuum limits. 
Finally, and equally importantly, one needs to know the relation of the coupling to 
a^'^ , the quantity of final interest, accurately to a few percent. This programme 
has been achieved for the pure gauge theory [H Ej. In full QCD calculations 
with Wilson fermions the amount of lattice data was barely enough to enable 
a reliable chiral and continuum extrapolation |21 Ej- Recent calculations with 
staggered fermions cover a wider range of lattice spacings and quark masses |3]. 
However, staggered fermions are not without their own problems. 

We determine A in the MS scheme from the force parameter tq jS] and the 
'boosted' coupling gn- The latter is obtained from the average plaquette. The 
advantage of this method is that both quantities are known to high precision. As 
in our previous work we shall use here non-perturbatively 0{a) improved 

Wilson (clover) fermions. Definitions of the action are standard (see, for example. 
Appendix D of [6 ). The lattice calculations will be done for n/ = 2 fiavours of 
dynamical quarks. In addition, we will update our quenched results. 

Since our first attempt [21 E] the amount of lattice data with dynamical quarks 
has greatly increased 7]. That is to say, at our previous couplings /? = 5.20, 
5.25 and 5.29 we have increased the statistics and done additional simulations 
at smaller quark masses. Furthermore, we have generated dynamical gauge field 
configurations at /? = 5.40 for three different quark masses. At each P value we 
now have data at three to four quark masses at our disposal, and the smallest 
lattice spacing that we have reached in our simulations is a ^ 0.07 fm. This allows 
us to improve on, and disentangle, the chiral and continuum extrapolations. In 
the quenched case the force parameter ro/a is now known up to /3 = 6.92 [Hj. 

The paper is organised as follows. In sectional we present a general discussion 
about the (3 function, including Fade approximations, and the running coupling 
constant. Also given are results in the MS scheme. In section IHl we set up the 
lattice formalism and discuss what coefficients are known. Various possibilities 
for converting to the MS scheme are given, which will indicate the magnitude 
of systematic errors. In section 0] results are given for TqA'^'^ for both quenched 
[uf = 0) and unquenched nj = 2 fermions. These results are then extrapolated 
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to Hf = 3 flavours of dynamical quarks in sectional This is done by matching 
the static force at the scale vq. In section we convert our results to physical 
units and, after matching to n/ = 5 flavours, compare them with other lattice 
determinations and to the experimental values. Finally, in section [3 we give our 
conclusions. 



2 The QCD Coupling and the (3 Function 



The 'running' of the QCD coupling constant as the scale changes is controlled by 
the P function, 

P^9s{M) (1) 



d\ogM 
with 

{gs) = -bogl - h9% - h^29l - big's - • • • , (2) 

renormalisation having introduced a scale M together with a scheme S. The first 
two coefficients are scheme independent and are given for the SU (3) colour gauge 
group as 

Integrating eq. ((T)) gives 

^ = F^{gs{M)) , (4) 



with 



1 1 bi 

+ 



P^iO bo^' bll 



(5) 

where A*^, the integration constant, is the fundamental scheme dependent QCD 
parameter. The integral in eq. (0) may be performed numerically or to low orders 
analytically. For example, to 3 loops we have 



M 

where 



A"" = h + ^Jhl - Abob^ 



(7) 
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Figure 1: versus /i/A*^^ for = (left picture) and Uf = 2 (right 

picture), using successively more and more coefficients of the /? function. 



and 



5 ^1 



hi - 26o&f 



bi 



" '1,2 ^~ 



(8) 



4bl^bl-4bobl 



Results are usually given in the MS scheme, with the scale M being replaced 
by /i, and thus 

A MS 

F«^(^?^(/i)) . 



In this scheme the next two (3 function coefficients are known ITIH ITT]: 
1 /2857 5033 325 

"2 



(9) 



blF - ^ 



nr H n"f I , 

18 ^ 54 ^ ' ' 



(4vr) 



149753 



6 



+ 3564C3 



+ 



162 

/ 50065 6472 



1078361 6508, , 



27 



(10) 



V 162 



C3 n^f + 



1093 



729 ^ 



The running coupling a^^{^) = g^{fi)/4:Tr is plotted in Fig. [T]for nj = 0, 2 by 
solving eq. (jSJ numerically, using only the first coefficient (1-loop), the first and 
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second coefficients (2-loop) etc. of the (5 function. The figure shows an apparently 
rapidly convergent series (cf the 3- to 4-loop result), certainly in the range we 
will be interested in, /i/A""^ ~ 20. The main difference between the = and 
Uf = 2 results is that af^|„^=2 rises more steeply as a function of /x/A"^, as 

bo\nf=2 < bo\nf=0- 

A knowledge of the (3 function to 4 loops is the exception rather than the rule. 
In many schemes it is known only to 3 loops. To improve the convergence of the 
P function, we may attempt to use a Pade approximation by writing eq. © as 



(^imM = } ,s I' ^ ' , (11) 



which on expanding is arranged to give the ffist three coefficients of eq. and 



estimates the next coefficient 6f as 



&f - ^ . (12) 
It is again possible to give an analytic result for using iif-yixy ^^"^ 



7V5 _ 



(13) 



At least for the MS scheme this appears to work reasonably well. Equation ()12p 
gives 63^^ ~ 3.22 x 10~^ and 1.67 x 10~^ for quenched and unquenched fermions, 
respectively, to be compared with the true values from eq. (fTIH) of 4.70 x 10~^ and 
2.73 X 10~^. In 12] we have shown a figure of the various Pade approximations to 
the /3 function. In Fig. |21we show the value of F"^{gjjs) at = 2 versus the 
/3 function coefficient number for both quenched and unquenched fermions. Also 
shown are the results using the [1/1] Pade approximations. It is seen that these 
numbers lie extremely close to the 4-loop (3 function results. As Pade approxi- 
mations give some estimation of the effect of higher order jS function coefficients, 
we shall thus prefer these later in our determination of the A parameter. 



3 Lattice Methods 

On the lattice we also have a coupling constant go{a) and corresponding (3 func- 
tion with coefficients bf^ and parameter A^, where 

aA^ = F^{go{a)). (14) 
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Figure 2: F"^(fi'37s) for = 2 versus (3 function coefficient number n. The 
Uf = values are filled circles, while the n/ = 2 values are filled squares. The 
[1/1] Fade approximations are given as open symbols. 



To evaluate , we need to know the 6f^s. They can be found by expanding 
gws as a power series in go as 

1 1 



9 ms 



+ 26o In a/i - tf^ + (26i In a/i - tf^)^^(a) 



+ [-260&1 In^ a/i + 2(6f ^ + hitf) Ina/x - t^] g^{a) + 
To have consistency between eqs. ^ and (fT^ we need 



and 



'^3 



"2 



26n In 



00^2 



bf' + 26f ^tf^ + h{t^y - 2bot^ 



(15) 



(16) 



;i7) 



where bf^ are the lattice /3 function coefficients, as in eq. So the transfor- 
mation between the two schemes is given by the (which define the transfor- 
mation), and the renormalisation group dictates how the scale running occurs 
(in this case the Ina/x terms). A knowledge of (the 1-loop) determines the 
relationship between the A parameters in the two schemes, while also knowing 
(the 2-loop) means that the 3-loop (3 function coefficient b^ can be found. 
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At present, what we know isOCSllIlliailHlinill 

tf^ = 0.4682013 - Uf [0.0066960 - 0.0050467 c,^ + 0.0298435 c^^ 

+ amg(-0.0272837 + 0.0223503 c,^ - 0.0070667 c^) + ^((amg))^] , 

(18) 

= 0.0556675 - n/[0.002600 + 0.000155 c,^ - 0.012834 
- 0.000474 - 0.000104 c^^ + O(amg)] . 

Here has been calculated including the arriq terms (m^ being the bare quark 
mass), while is known only for arriq = 0, and tif^ is unknown, which means 
that from eq. ()17|1 62''^ is known but not . For general Csw the connection 
between and is only defined up to terms of 0(a), but on the improvement 
trajectory Cgw = ^ + 0{gQ) it is possible to arrange it to be O(a^) if the anig terms 
are included in the tf^s. 



Thus, the conversion from the lattice coupling to the MS coupling (eqs. (fT3jl 
and ()18|) ) can also be written with mass independent tf^s, if we redefine g^ by 
replacing it by g^, where 

~gl = glil + b,am,) , b, = bfu^gl + 0{gt) . (19) 

So, putting Csw = 1 + 0((7q) into eqs. (fTH|l and (fTK|) means that t'{^ is replaced by 
tf^ — Hf amqb^g \ which gives bf^ = 0.01200. This value agrees with the number 
reported in [TH] . 

Thus, in this mass independent scheme (ie a scheme where the renormalisation 
conditions are imposed for zero quark mass) there appears to be little difference 
in extrapolating to the chiral limit using constant (3 = Q/gQ, rather than constant 
P = G/go^. So, rather than using eq. (fTHj) at finite arUg, we shall first extrap- 
olate our plaquette and ro/a data to the chiral limit and then determine A*^'^. 
Before attempting this, we shall discuss some improvements to help improve the 
convergence of the power series (jl5|) . 

As it is well known that lattice perturbative expansions are poorly convergent, 
we have used a 'boosted' coupling constant 

95 - # (20) 

H 

to help the series (fT^. or equivalently Q for P^'^{go), converge faster. Here 
P = Uq = (Trt/°)/3 is the average plaquette. In perturbation theory we write 

\ = \-Pi-P2g'o + 0{g'o) (21) 

with [1111201 



3 ' (22) 

) 



P2 = 0.0339110 - n/(0.001846 - 0.0000539 c,^ + 0.001590 
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for massless clover fermions. 

To improve the convergence of the series further, we re-express it in terms of 
the tadpole improved coefficient 

'^sw ~ CswUq ■ (23) 

Changing tf^ to t° first replaces tf^ by tf^ — Pi, and secondly using simply 
replaces every Cgw by in tf*^, but the change in ^3 is more complicated as the 
coefficients of change in ■ 

This gives for tf = tf{c°^) in the chiral limit 

t° = 0.1348680 - n/[0.0066960 - 0.0050467c°^ + 0.0298435 {c°J^] 
t° = 0.0217565 - n/ [0.000753 + 0.001053 - 0.000498 {c°J'^ (24) 
- 0.00047 4(c°^)3 - 0.000104 (c°J^] . 

As we have here a 2- loop result, we can see how well tadpole improvement im- 
proves the series convergence. The coefficient of n/ in t° is considerably smaller 
than the corresponding coefficient in . For example, using the values at 
j3 = 5.40 given in the next section, we find that the magnitude of the coeffi- 
cient is reduced by two orders of magnitude (from ~ —0.0438 to ~ 0.0003). 
What this tadpole improvement represents is taking a path from = to 
= (?□, keeping fixed. Later we shall consider other trajectories from 
to (yfp. If we had all orders of the theory, the result would depend only on the 
end point. But with a finite series the trajectory will matter. This will help us 
estimate systematic errors from unknown higher order terms. 
Thus, in conclusion we have 

aA° = F°(^o(a)), (25) 



together with the conversion formula 



A MS 

F^%9M^^)) , (26) 



with 



^ + 2bQln afi - t° + {2bilnai^-t°) glia) + .. . (27) 



A MS 

i° = 26oln^ (28) 



and 

6° = bf^ + bit° bot° . (29) 



We shall now discuss various strategies to determine A 



MS 
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3.1 Method I 

This method was used in our previous papers |21IS1I2I], with the difference that 
now we first extrapolate to the chiral hmit. For each j3 value we first compute 
tf from eq. (j^ . Then from eq. ^Tl\ we convert ga to gjjs at some appropriate 
scale /i*, and using the force scale tq, we calculate tqA"^ from eq. 



roP^'" = ro^i.F'^''{g^s{^^,)). (30) 

Finally, we extrapolate to the continuum limit, a 0. Note that t° will depend 
on the coupling because does. 

We must determine the scale /x*. A good choice to help eq. (jTTj) converge 
rapidly is to take the 0(1) coefficient to vanish, which is achieved by choosing 

jjL^ = -exp TTT- . (31) 



a \2bo 

Thus, we used 



+ {hi-t2]9M+0{g',) (32) 



to find (7|j-g-(/i=K), which was then substituted into eq. (jHUjl . 



3.2 Method II 

Alternatively, we can first determine 6° from eq. ()29p and then determine roA^ 
via eq. After computing this, we convert to tqA*^^ using 

roA^ = roA°exp(^^^, (33) 

and then take the continuum limit. Again, note that 63 ^i^^ depend on the 
coupling, because does. 

This method is equivalent to choosing a scale /i=, as in method I, such that 
gMs{fJ^=) = gn{a). In this case all the coefficient terms of eq. (jTfjl vanish. The 
scale that achieves this is 

Indeed, substituting /i= into eq. then gives eq. again. The scale //= is 
close to yU*, as can be seen by expanding eq. (fH^ to 3 loops. From eq. (jHI) we 
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have 



1 ft^ 



1 + ^g' 



-Pa 



1 + 



-PS 



Pa SMS 2^ 



1 1 ^I'^l ~ ^0*^2 2 I 
1 + • • • 



26g 



> At* , 

for the couphngs used here. 



(35) 



3.3 Method III 



Another possibihty, and theoretically the most sound, is to vary along the 
improvement path as g"^ increases. This will give genuinely constant 13 function 
coefficients (ie independent of the coupling). As the 1-loop expansion for is 
known along this path, 

= 1 + Co + . . . , (36) 
with Cq = Co — |pi and Cq = 0.2659(1) 22^, then expanding eq. (pij) gives 



^0^2 I =1 



boCo 



-0.0008241 



(37) 



This result may also be derived from eq. ()27|) by first setting a = /i ^ (for 
simplicity) and then taking fid/dfi of this equation. This leads to 



9ms 



9n dc°^ dgn 



(38) 



which upon expanding out also gives eq. 

So, having determined 62 ^q. (j29p . the method is as for method II: first 
determine roA° using eq. and then convert to roA"-^ using eq. (j2Sl)- 



3.4 Methods IIP and HIP 

To further improve our calculations, and to reduce the systematic error, we con- 
sider here the effect of Pade improving the /3 function, as given in eqs. (|TT|l and 
(IT!?jl . We restrict ourselves to methods II and III, and we call the Pade improved 
results IIP and HIP, respectively. 
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u.ovoyj t yi uo i 
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U . U U / (J / O/f J 


6.07 


6.033(17) 


0.601099(18) 


0.5696(16) 


0.5746(16) 


0.5692(16) 


6.20 


7.380(26) 


0.613633(02) 


0.5861(21) 


0.5907(21) 


0.5855(21) 


6.40 


9.740(50) 


0.630633(04) 


0.5976(31) 


0.6018(31) 


0.5970(31) 


6.57 


12.18(10) 


0.643524(15) 


0.6029(48) 


0.6067(48) 


0.6022(48) 


6.69 


14.20(12) 


0.651936(15) 


0.6055(50) 


0.6091(51) 


0.6049(50) 


6.81 


16.54(12) 


0.659877(13) 


0.6080(46) 


0.6113(46) 


0.6073(46) 


6.92 


19.13(15) 


0.666721(12) 


0.6145(47) 


0.6111(41) 


0.6139(41) 


oo 


oo 


1 


0.6152(21) 


0.6189(21) 


0.6145(20) 



Table 1: The quenched roA""^ values for methods 1, II and IIP (ie using the Pade 
improved (3 function /5[i/i]) together with the force parameter ro/a |H] (the number 
dX (3 = 6.0 is from the interpolation formula given there) and the plaquette P. 
The continuum extrapolated values together with the statistical errors are given 
in the bottom row. Numbers in italics are not used in the fits. 
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Figure 3: The quenched tqA"^ points versus (a/ro)^, together with a linear ex- 
trapolation to the continuum limit for method IIP. The filled circles are used for 
the extrapolation. The star represents the extrapolated value. 
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4 Results 



4.1 Quenched Results 

In the quenched case [uf = 0) we do not have any of the additional chiral hmit 
extrapolation complications alluded to in the previous section, or a Cg^ term. 
This means that there is no difference between method II and method III, so 
the procedure is straightforward. In Table ^ we give the parameters used. For 
ro we use, for consistency, exclusively the values given in [Hj, which includes 
previous results from The one exception is P = 6.0, where we have used the 
interpolation formula [HI for rg/a. Our plaquette values are determined at their 
given P values. 

In Tabled we also give the results for tqA"^ from methods I, II and IIP. We 
first see that the results for tqA*^^ are almost indistinguishable between methods 
I, II and IIP. Method IIP lies just below method I (and indeed is almost identical 
with it). 

We now consider the continuum limit of our results. In Fig. 01 we plot the 
results for tqA""^ against (a/ro)^ for method IIP. The differences between the 
results of the various methods are small. As one expects that Pade improvement 
gives a better answer, we shall concentrate on IIP. The smallest a value is not 
included in the fit, as it appears to deviate a little, but including it would not 
have changed the extrapolated value much. We also have not included /3 = 6.0 in 
the fit, as To/a is only known from an interpolation formula. But as can be seen 
from the figure, including it has no effect on the result. Also, the two coarsest 
a values have not been included in the fit, as they show significant non- linear 
effects in a^. These two points are not shown in the plot, as they lie far to the 
right. Figure El clearly shows a linear extrapolation over a wide range of lattice 
spacings, ~ 2 - 6.5 GeV, giving a value for method IIP of 

roAr = roA^L^=o = 0-614(2) (5). (39) 

Here the first error is statistical, and the second systematic error is estimated by 
the spread in the results between methods I, II and IIP. That the systematic error 
is small is an indication of the convergence of results from the different methods. 
The result (jHUj) agrees with our earlier value P|. 

4.2 Unquenched n/ = 2 Results 

We now turn to unquenched Uf = 2 fermions. In Table El we show the (3, k and 
Csw parameters used in the simulations, together with the measured vq/q and 
plaquette P values. As discussed in section ^ we shall first determine tqA^''^ in 
the chiral limit and then perform the continuum extrapolation. We must thus 
first find the zero quark mass results from Table El We shall make a chiral 
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K 


V 




ro/a 


F 


Group 


5.20 


0.1342 


16^ X 32 


2.0171 


4.077(70) 


0.528994(58) 


QCDSF 


5.20 


0.1350 


16 X 32 


2.0171 


4.754(45) 


0.533670(40) 


UKQGU 


5.20 


0.1355 


16 X 32 


2.0171 


5.041(53) 


0.536250(30) 


UKQGD 


5.20 


0.13565 


16 X 32 


2.0171 


5.250(75) 


0.537070(100) 


UKQGD 
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5.25 


0.1346 


16^ X 32 


1.9603 


4.737(50) 


0.538770(41) 


QGDSF 


5.25 


0.1352 


16^ X 32 


1.9603 


5.138(55) 


0.541150(30) 


UKQGD 


5.25 


0.13575 


243 X 48 


1.9603 


5.532(40) 


0.543135(15) 


QGDSF 


5.29 


0.1340 


16^ X 32 


1.9192 


4.813(82) 


0.542400(50) 


UKQGD 


5.29 


0.1350 


16^ X 32 


1.9192 


5.227(75) 


0.545520(29) 


QGDSF 


5.29 


0.1355 


243 X 48 


1.9192 


5.566(64) 


0.547094(23) 


QGDSF 


5.29 


0.1359 


243 X 48 


1.9192 


5.880(100) 


0.548286(57) 


QGDSF 


5.40 


0.1350 


243 X 48 


1.8228 


6.092(67) 


0.559000(19) 


QGDSF 


5.40 


0.1356 


243 X 48 


1.8228 


6.381(53) 


0.560246(10) 


QGDSF 


5.40 


0.1361 


243 X 48 


1.8228 


6.714(64) 


0.561281(08) 


QGDSF 



Table 2: The unquenched f3, k and Csw values and the volume V, together with 
the measured force parameter ro/a and plaquette P. The collaboration that gen- 
erated the configurations is given in the last column. The results for f3 = 5.29, 
K = 0.1359 are preliminary. We have reanalysed our ro/a values, taking auto- 
correlations properly into account, which gave larger error bars than previously 
reported ^24j. 







ro/a 


P 


5.20 


0.136008(15) 


5.455(96) 


0.538608(49) 


5.25 


0.136250(07) 


5.885(79) 


0.544780(89) 


5.29 


0.136410(09) 


6.254(99) 


0.549877(109) 


5.40 


0.136690(22) 


7.390(26) 


0.562499(46) 



Table 3: The critical values for n (ie k^) and the chiral limit values for ro/a and 
P for the four /? values used here. 
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Figure 4: The plaquette P (filled symbols) plotted against the bare quark mass 
arUg for (3 = 5.20 (lower curve) until P = 5.40 (upper curve). The fits use eq. (jHJ, 
giving the extrapolated values in the chiral limit (open symbols). 
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Figure 5: The force parameter ro/a plotted against arriq. The same notation as 
in Fig. m is used. 
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/5 


TqA^ I 


roA^ II 


roA^ IIP 


roA^ III 


roA^ HIP 


5.20 
5.25 
5.29 
5.40 


0.5183(91) 
0.5210(71) 
0.5372(85) 
0.5577(198) 


0.5304(94) 
0.5415(73) 
0.5482(87) 
0.5676(201) 


0.4913(87) 
0.5040(68) 
0.5120(81) 
0.5343(189) 


0.6459(114) 
0.6450(87) 
0.6433(102) 
0.6431(228) 


0.6173(109) 
0.6174(83) 
0.6165(98) 
0.6182(219) 


oo 


0.6012(346) 


0.6085(352) 


0.5819(329) 


0.6376(412) 


0.6170(395) 



Table 4: The values for tqA^"^ for methods I, II, IIP, III, HIP described in sectional 
for the four P values used here. 

extrapolation in arriq, defined here by 

We estimate Kc from partially quenched pion data. The results have been given 
in [7j and are tabulated in the second column of Table IHl 

In Fig. m we show the results for the plaquette. For each j3 value the data 
appear to be rather linear in the quark mass arriq, in particularly for the higher /3 
values. This suggests that a quadratic fit ansatz is sufficient to obtain the value 
of P in the chiral limit. We thus use 

P = do + diairiq + d2{amgY . (41) 

Except for (3 = 5.20, it does not much matter whether a linear or quadratic fit is 
used. In Fig. we show the results for rg/a. The data is less linear in arriq, and 
also less smooth, so we used the renormalisation group inspired global fit ansatz 

In ^ = A,{i3) + A2{f3)amq + A^{l3){amqf , (42) 

where ^i(/9) is a linear polynomial in and A2{I3), A2,{(3) are quadratic polyno- 
mials in 13. This ansatz was also used in 0. The results of the fits in the chiral 
limit are given in Table 121 

In Table |3] we give our results for roA""^ for methods I, II, IIP, III and HIP. 
Again, as the results for method I are very similar to method II, we shall not 
discuss method I further here. In Fig.Elwe plot tqA.^'^ against {a/r^Y for methods 
IIP and HIP, together with a linear extrapolation to the continuum limit. Though 
we cannot reach such small a values as for the quenched case, the tqK'^'^ data do 
seem to lie on straight lines. We find a linear behaviour at least over the region 
~ 2 - 3 GeV. This seems to be well inside the linear region of Fig. El 

For methods IIP (and II) the results lie roughly parallel to the quenched 
results, while for methods HIP (and III) they are fiatter and higher. However, 
in the continuum limit they agree within error bars. Ideally, the result should 
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Figure 6: The unquenched tqA"^ points (filled circles) versus (a/ro)^, together 
with a linear extrapolation to the continuum limit for methods IIP and HIP. 
Stars represent the extrapolated values. 

not depend on the choice of trajectory. The way this should work, as mentioned 
before, is that although the coefficients tf will be different depending on the path 
one might choose, the sum 



should not. However, at the order to which we have the series this is not yet so. 
The difference between methods II and III is that we have replaced by its 
1-loop expansion. Returning to Fig. IHl the fact that the results from methods II, 
IIP are almost parallel to the quenched results suggests that in methods II, IIP 
the O(a^) effects come from the same source as in the quenched case, which must 
be the gluon action. For methods III, HIP the slope is much smaller so there 
must have been a fortuitous cancellation between effects from the gluon and 
fermion terms. 

One expects that Pade improvement gives a better answer, so the P results 
are more trustworthy. Previous experience suggests that the procedure in IIP of 
using tadpole improved Csw works fairly well. For example, Kc in |23 and the 
renormalisation constant Z for V2b in agree within a few percent with the 
non-perturbative values. However, method HIP is a more consistent approach. 
Furthermore, the results from method HIP appear to be insensitive to the par- 
ticular form of the continuum extrapolation. We therefore take these numbers as 



1 



t°-t°^?n(a) + ... 



(43) 
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our best estimate. 

From the linear extrapolation of method HIP to the continuum limit we thus 
quote 

roAr = = 0.617(40)(21) , (44) 

where the first error is statistical and the second systematic. The latter error is 
estimated by the spread in the results between method III and HIP. Compared 
to our previous result |2j, the value (jl^ has increased by ~ 10%, but still lies 
within the error bars. 



5 Extrapolation to n/ = 3 Flavours 

At high energy scales we can see that A'*^^ makes some fairly large jumps as we 
pass through the heavy quark mass thresholds and change the effective number 
of flavours. From 120] we can see that the reason for these large jumps is the 
fact that TTiq/A*^^ is large. We want to argue here that the situation with light 
quarks, < A^'''^, is rather different, and that in this case we do not expect to 
see any dramatic dependence of A*^'^ on n^. 

We will determine the Uf = 3 flavour A parameter from matching the static 
force at the scale tq. 



5.1 One- loop Matching 

To make clear what is involved in matching, we will go through the 1-loop cal- 
culation in some detail. 

At the 1-loop level the static potential between fundamental charges is given 

by 



3 Airr 



31 

22 ( ln/ir + 7£; + — 



4 / 5 
-Uf I In /ir + 7£ + - 



(45) 



+ 



for massless sea quarks (see, for example, [ZZ])- We can work out the force /(r) 
at distance r by differentiating this to give 



4vrrV(r) = ^^L(/.)|l + ^ 



22 In fir + -^E 



35 
66 



4 A 1 



(46) 



+ 
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If we now change the flavour number from 2 to 0, or from 2 to 3, while keeping 
the force at distance r constant, we get 

(33-6)ln-^ = 2^1nArr + 7,.-- 
We can eliminate r from these equations, leaving us with the simple equation 

Ap ~ VAp ) ' 
which can be used to estimate A3^^ from the Uf = and Uf = 2 results. 




5.2 Higher Loops 

To repeat this matching calculation with more loops, we follow jH] and define a 
force-scale coupling g^^ by 



(49) 



From eq. ()46|) we can read off 

1 



(4vr) 



35\ 4 / 1 



(50) 



We can find tg' by calculating the force from the 2- loop expression of V{r) re- 
ported in |28l^: 



-l-iq 
''2 



(An) 



1107 , 229 . 9 4 

204 vr^ + -vr^ - 66 Cs 

2 '34^^ 



Uf 



^ + 767.. + + 52 Cs) + ^nj (12 - vr^) 



(51) 



which gives us enough information to calculate the 3-loop j3 function for gqg{r) (cf 
eq. (fT7|)). There would be complications in going to the next order, because it is 
known that terms of the type Ina^ will enter the series for the potential jHUj . 

We are now ready to see how A^-'^ depends on flavour number, if we make the 
value of /(r) independent of (the number of massless quark flavours) at some 
particular r value. Implicitly, we assume r ^ vq. If /(r) is independent of Uf, 
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Figure 7: The ratio A3^'^/A"'^ against the ratio A^^/Aq^ from 1-loop (dotted 
hne), 2-loop (dashed hne) and 3- loop (sohd hne) matching. 

then gqq{r) is independent of n/ too. We can compare the qq scheme As by using 

rA^^ = F'%,,(r),n; = 0), 

rAf = F^^{g,,{r),nf = 2), (52) 
rAf = F^^g,,{r),nf = 3). 

We can take ratios of these equations to cancel r and find equations for A ratios. 
These qq scheme A ratios can then be converted into MS by using tf^ from 
eq. (jSUj) . This gives us a way of making a parametric plot of A ratios by varying 

and calculating all three As from g^^. In Fig. [7|we show the plot. 

The results clearly have to be treated with some caution, because r^A is a 
fairly large number. So it is not clear how much we can learn from perturbative 
results at the scale vq. It is therefore quite surprising that the different orders 
of perturbation theory agree so well in Fig. [7| Furthermore, we have assumed in 
this section that rorus <^ 1, so that the strange quark can reasonably be treated 
as massless. Both these difficulties could be decreased by using a smaller distance 
(and thus a smaller value for r^/(r)) to set our scale. 

5.3 Result for n/ = 3 

From our quenched and unquenched nj = 2 results, ()39p and ()44|) . we obtain 
A"'^/Ag^'^ = 1.005. If we insert this number into the 3-loop matching curve 
shown in Fig. [7| we find A"^/A2^^ = 0.999. From this ratio and eq. (j44j) we then 
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obtain for Uf = 3 quark flavours 

roAr = ^oA"^L^,=3 = 0.616(29)(19) . (53) 

We liave not attempted to estimate the systematic error induced by the matching 
procedure. 

6 Comparison with Phenomenology 

In this section we shall make a comparison with other lattice and phenomeno- 
logical results. For this we first need to set the force scale in terms of a physical 
unit. 

To fix the scale ro in physical units, we extrapolate recent dimensionless nu- 
cleon masses m^ro found by the CP-PACS, JLQCD and QCDSF-UKQCD col- 
laborations jointly to the physical pion mass following [21]. This gives the value 
vq = 0.467 fm with an estimated error of 7%. We will use this number throughout 
this paper. A similar result for tq was quoted in [3T] . 

For the quenched case we then obtain 

A^ = 259(1) (20) MeV, (54) 

and for the unquenched case we find 

Af^ = 261(17)(26)MeV, (55) 
Af^ = 260(12) (26) MeV. (56) 

The systematic errors quoted here include the uncertainty in setting the scale. 
Note that previously we had assumed vq = 0.5 fm. 

In Fig.lSlwe show our results for A"^ together with recent experimental values 
from and J^]. It appears that the lattice results extrapolate smoothly to the 
experimental values at n/ = 4 [3^ and nj = 5 33^. However, our rij = 3 re- 
sult lies two standard deviations below the corresponding phenomenological value 
(open triangle). (The reader should be aware that the sometimes called experi- 
mental numbers imply a good deal of modelling and, thus, should be regarded as 
phenomenological numbers.) 

In order to compare as from various experiments and theory, it must be 
evolved to a common scale. For convenience this is taken to be the mass of the 
Z boson, mz- Having computed A*^® for n/ = 3 fiavours, we may use the 4-loop 
expansion of as and the 3-loop matching condition at the quark thresholds PSIIH^ 
to determine a^^^^{mz). We take the charm and bottom thresholds to be at 1.5 
and 4.5 GeV, respectively. Furthermore, we choose the charm and bottom quark 
masses to be m'^^{m^ = 1.5 GeV and mff '^lmb) = 4.5 GeV, respectively. Varying 
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Figure 8: Values of A^^'^ versus number of quark flavours Uf. The filled circles are 
our Hf = 0, 2 results, and the open circle is our extrapolated value. The inner 
error bars give the statistical errors, while the outer error bars give the total errors. 
The square is from a 3- loop analysis of the non-singlet structure functions |32j . 
The triangles are taken from The open triangles are evaluated using the 

4-loop expansion of as and 3-loop matching at the quark thresholds. The entries 
at n/ = 3 and 4 have been displaced horizontally. 



the charm and bottom quark masses within reasonable limits has a neglible effect 
on the final result. We then obtain 

<=5M = 0.112(1)(2). (57) 

This is to be compared with the world average value [32] al''^{mz) = 0.1182(27). 

In Fig. iniwe compare our result for af^{mz) with other lattice results and 
experiment. We find agreement with previous lattice calculations using Wilson 
fermions. It occurs that the Wilson results lie systematically below the mean 
experimental value. On the other hand, calculations using staggered fermions 
(albeit from the same group) show a better agreement with experiment. Our 
result for rQA2^ agrees also with that of the ALPHA collaboration |3U], which 
does not quote a number for a^'^^mz)- Our result for a^^^{mz) lies two standard 
deviations below the phenomenological value. 
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Figure 9: Comparison of ag'^lmz) from this work (solid circle) with other lat- 
tice results ISaEElEZliaiMillMlini (from top to bottom). The circles are 
from Wilson fermions and the squares from staggered fermions. The dashed line 
indicates the mean experimental value 

7 Conclusions 

Due to substantial improvements of the performance of our hybrid Monte Carlo 
algorithm |1T], we were able to extend our dynamical simulations to smaller 
quark masses and to larger values of (3. Our smallest lattice spacing now is 
a ~ 0.07 fm. This enabled us to perform a chiral and continuum extrapolation 
of the lattice data. Because the calculation involves a perturbative conversion 
from the lattice coupling constant to the (mass independent) MS constant, it 
was important to first extrapolate the lattice data to the chiral limit. We have 
discussed basically two approaches of converting the lattice coupling constant to 
the MS one. They differed mainly in how the non-perturbative improvement 
(clover) term was incorporated in the perturbative expansion. It was reassuring 
to see that both methods led to the same result in the continuum limit. This 
indicates once more that a reliable extrapolation to the continuum limit is very 
important. 

We could also improve on our quenched result, because data at smaller lattice 
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spacings became available. 

There are several sources of systematic error in our calculation. The main 
error comes from setting the scale, followed by the continuum extrapolation. As 
better dynamical data become available, the uncertainty in setting the scale will 
be gradually reduced. Simulations at smaller lattice spacings will become possible 
with the next generation of computers, which should facilitate the extrapolation 
to the continuum limit. 
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